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We investigate the thermally activated magnetization 
switching of small ferromagnetic particles driven by an ex- 
ternal magnetic field. For low uniaxial anisotropy the spins 
can be expected to rotate coherently, while for sufficient large 
anisotropy they should behave Ising-like, i. e. , the switching 
should then be due to nucleation. We study this crossover 
from coherent rotation to nucleation for the classical three- 
dimensional Heisenberg model with a finite anisotropy. The 
crossover is influenced by the size of the particle, the strength 
of the driving magnetic field, and the anisotropy. We discuss 
the relevant energy barriers which have to be overcome during 
the switching, and find theoretical arguments which yield the 
energetically favorable reversal mechanisms for given values 
of the quantities above. The results are confirmed by Monte 
Carlo simulations of Heisenberg and Ising models. 
Pacs: 75.10.Hk, 75.40.Mg, 64.60.Qb 



I. INTRODUCTION 

The size of magnetic particles plays a crucial role for 
the density of information storage in magnetic record- 
ing media. Sufficiently small particles become single- 
domain particles, which improves their quality for mag- 
netic recording. On the other hand, when the particles 
are too small they become superparamagnetic and no 
information can be stored (see e. g. Q for a review). 
Hence, much effort has recently been focused on the un- 
derstanding of small magnetic particles, especially since 
recent experimental techniques allow for the investigation 
of isolated single-domain particles [||-^| . 

In this paper focus is on the reversal of ferromagnetic 
particles of finite size. We investigate the influence of the 
size and anisotropy of the particle on the possible reversal 
mechanisms, two extreme cases of which are coherent ro- 
tation and nucleation. The latter mechanism has been a 
subject of common interest in recent years |3|-[n], stud- 
ied mainly theoretically in Ising models, which can be 
interpreted as a classical Heisenberg model in the limit 
of infinite anisotropy. It is the aim of this paper to study 
the crossover from magnetization reversal due to nucle- 
ation for high anisotropy to coherent rotation 12 ll| for 
lower anisotropy. 
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Throughout the paper we will consider a finite, spher- 
ical three- dimensional system of magnetic moments. 
These magnetic moments may represent atomic spins or 
block spins following from a coarse graining of the phys- 
ical lattice |l4j]. Our system is defined by a classical 
Heisenberg Hamiltonian, 



where the Si are three-dimensional vectors of unit length. 
The first sum which represents the exchange of the spins, 
is over nearest neighbors with the exchange coupling 
constant J. The second sum represents an uniaxial 
anisotropy which favors the z-axis as easy axis of the 
system (anisotropy constant d > 0). The last sum is the 
coupling of the spins to an applied magnetic field, where 
B is the strength of the field times the absolute value of 
the magnetic moment of the spin. We neglect dipolar in- 
teraction. Although in principle it is possible to consider 
dipolar interaction in a Monte Carlo simulation [H this 
needs much more computational effort due to the long 
range of the dipolar interaction and, hence, exceeds cur- 
rent computer capacities . Therefore, the validity of our 
results is restricted to particles which are small enough 
to be single-domain in the remanent state Jig] . 

In the following, we will investigate the thermally ac- 
tivated reversal of a particle which is destabilized by an 
magnetic field pointing in a direction antiparallcl to the 
initial magnetization which is parallel to the easy axis of 
the system. Due to the finite temperature and the mag- 
netic field, after some time the particle will reverse its 
magnetization, i. e. the z-component of the magnetiza- 
tion will change its sign. 

In section |ll|, we determine the energy barriers which 
have to be overcome by thermal fluctuations for the two 
cases of coherent rotation and nucleation within a clas- 
sical theory. By a comparison of the energy barriers, we 
derive where the crossover from one mechanism to the 
other occurs. 



In section III we compare our theoretical considera- 



tions with numerical results from Monte Carlo simula- 
tions of Heisenberg and Ising models, and we relate the 
lifetime of the metastable state to the theoretical energy 
barriers for the different reversal mechanisms. 



II. THEORY 

A. Coherent rotation 

We give here a brief summary of the results of the 
theory of Neel |Q and Brown Q, since we need this 
concept for the further progress of our theoretical con- 
siderations. 

Let us consider a spherical, homogeneously magnetized 
particle of radius R. The simplest theoretical description 
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for the reversal of such a particle is to assume that the 
reversal mechanism is coherent rotation, i. e. a uniform 
rotation of all spins of the particle. This reversal process 
can be described by an angle of rotation 9 between the 
easy axis of the system - which in our case will be an- 
tiparallel to the direction of the magnetic field - and the 
magnetization of the particle. The increase of the energy 
during the reversal is then 



AE = -dV cos 2 9 - BMV cos 9. 



(2) 



Since this equation should be comparable to Eq. [IJ the 
anisotropy constant d is an anisotropy energy per unit cell 
(spin), and V — (4/3)7ri? 3 is the volume of the particle 
as a number of unit cells. B is - as before - the absolute 
value of the applied field times the magnetic moment of 
a unit cell and, hence, M the spontaneous magnetization 
per magnetic moment. The energy barrier which has to 
be overcome is due to the anisotropy of the system. It is 
the maximum of AE with respect to 9: 
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The corresponding lifetime of the metastable state is then 
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for temperature T <C AE cr . The two equations above 
are physically relevant only for 



d > MB/2. 



(5) 



Otherwise there is no (positive) energy barrier, and hence 
the reversal is spontaneous without the need for thermal 
activation. This is the region of nonthermal reversal. 



B. Nucleation 




FIG. 1. The nucleation of a droplet at the boundary of a 
spherical particle. 

Then the surface of the domain wall is F = 2nrh and 
the volume of the droplet (the shaded region in Fig. |l|) 
is V d = nH 2 (3R - H)/3 + Trh 2 (3r - h)/3. The energy 
increase during the reversal of the particle is 



AE = 2aF - 2MBV d , 



(6) 



where a is the energy density of the domain wall. Fur- 
thermore, it is h = r — \J r 2 — 2HR + H 2 , since the quan- 
tities in Fig.[l] are not independent. Hence the energy in- 
crease can be expressed in terms of H , which is a measure 
for the penetration depth of the domain wall and its cur- 
vature r. These two quantities define the geometry of the 
droplet. Next we determine that curvature r mm which 
minimizes the energy increase by the condition 
yielding the physically relevant solution 
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which is also the radius of a critical droplet for classical 
nucleation in a bulk material. 

The energy barrier AE W which has to be overcome 
during the reversal is the maximum of the energy increase 
with respect to H. From the condition = 0, the 

physically relevant solution 



For a system with a sufficient large anisotropy, it might 
be energetically favorable to divide into parts with oppo- 
site directions of magnetization parallel to the easy axis 
in order to minimize the anisotropy energy barrier. This 
kind of reversal mechanism is called nucleation |ll| (see 
Q for a recent review). The simplest case of a reversal 
process driven by nucleation for a system of finite size 
is the growth of one single droplet starting at one point 
of the boundary of the system (see also llTJ for a corre- 
sponding calculation for a two-dimensional system) . Due 
to the growth of the droplet a domain wall will cross the 
system, and the energy barrier which has to be overcome 
is caused by the domain- wall energy. We assume that the 
domain wall will have a curvature defined by a radius r 
(see Fig. |). 



= R 1- 
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follows with x = MBR/a. Inserting the two conditions 
above into the formula for the energy increase yields the 
energy barrier for nucleation 



AE n = 
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This expression has two important limits. The first is the 
limit of infinite system size, 



lim AE n = 
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where we obtain half of the energy barrier of the classical 
nucleation theory for a bulk system. The reduction by 
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a factor of \ is due to the fact that for open boundary 
conditions only one-half of a critical droplet has to enter 
the system from the boundary. 

The other interesting limit is that of small magnetic 
fields, where Eq.|] can be expanded with respect to x = 
MBR/a resulting in 



AE n « 2ttR 2 <7 
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This means that for a small field B the energy barrier 
of a nucleation process is the energy of a flat domain 
wall in the center of the particle plus corrections which 
start linearly in B. In contrast to the Neel-Brown theory, 
here, for vanishing magnetic field, the energy barrier is 
proportional to the cross-sectional area of the particle 
rather than its volume (see also the work of Braun jl7| ) , 
which consequently reduces the coercivity of the particle. 



C. Comparison of coherent rotation and nucleation 

Comparing the two energy barriers for coherent rota- 
tion (Eq.[j|) and nucleation (Eq.^) , we can evaluate which 
reversal mechanism has the lower activation energy for a 
given set of values of R, B, a, and d. A corresponding 
diagram is shown in Fig. || for a system of radius R = 4 
spins, where we set M — 1/spin and a = J/spin. 
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FIG. 2. Diagram showing the regions of different reversal 
mechanisms for a particle of size R = 4 spins, with a — J/spin 
and M = 1/spin. The points are results from Monte Carlo 
simulations (see Section ill). 



For large anisotropy the reversal is dominated by nu- 
cleation - the particle behaves like an Ising system. The 
crossover line d c (B,a,R) which separates the region of 
reversal by nucleation from the region of coherent rota- 
tion can be determined by the condition that here the 
energy barriers for a nucleation process and for coherent 
rotation are equal. This condition results in 
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For vanishing magnetic field, the formula above has a 
finite limit. 



lim (d c ) = !^r. 
b^o 2R 



(13) 



For large particle size, Eq. [12| has the simple asymptotic 
form 
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which, interestingly, is also the limit for nonthermal re- 
versal (see Eq. H) . This means that for increasing particle 
size the region where a thermally activated reversal by 
coherent rotation occurs vanishes. For an infinite parti- 
cle size the reversal is always either nonthermal or it is 
driven by nucleation - depending only on the ratio of the 
magnetic field to the anisotropy. 

All the considerations above can be expected to be rel- 
evant only for sufficiently low temperatures. For higher 
temperatures the situation is more complicated. That is, 
for the nucleation regime, a crossover from single-droplet 
to multi-droplet nucleation - with different energy bar- 
riers) - is discussed in the literature (see |T§|-|2"0|] and 
references therein). 

Apart from that, in a system with a given finite 
anisotropy the domain walls may be extended to a cer- 
tain domain wall width £. For large anisotropy £ becomes 
as small as one lattice constant. This is the Ising case, 
where the domain-wall energy density is a ~ J/spin. 
For smaller anisotropy the domain walls become more 
extended, and a decreases. Hence the crossover from 
nucleation to coherent rotation may be softened by the 
occurrence of extended domain walls. In this sense, pure 
coherent rotation could also be interpreted as a domain- 
wall-driven reversal, where the width of the domain wall 
is larger than the particle size. Obviously, our theoreti- 
cal considerations discuss only two extreme cases. How 
realistic they are has to be tested numerically. 



III. MONTE CARLO SIMULATION 
A. Method 

Due to the many degrees of freedom of a spin system, 
numerical methods have to be used for a detailed micro- 
scopic description of the system. Since we are especially 
interested in the thermal properties of the system, we use 
Monte Carlo methods pi ] for the simulation of the mag- 
netic particle. Although a direct mapping of the time 
scale of a Monte Carlo simulation on experimental time 
scales is not possible, this method provides information 
on the dynamical behavior of the system since it solves 
the master equation for the irreversible behavior of the 
system p2[ . 

We consider spins on a simple cubic lattice of size L x 
L x L, and simulate spherical particles with radius R = 
L/2 and open boundary conditions on this lattice. 
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One single spin flip of our Monte Carlo procedure con- 
sists of three parts. First, a spin is chosen randomly and 
a trial step is made (the role of which we will discuss 
below). Second, the change of the energy of the sys- 
tem is computed according Eq. [lj Third, the trial step 
is accepted with the probability from the heat-bath al- 
gorithm. Let us call one sweep through the lattice and 
performing the procedure explained above once per spin 
one Monte Carlo step (MCS). 

Since we are interested in different reversal mecha- 
nisms, we designed a special algorithm which can sim- 
ulate all of them efficiently. We use three different kinds 
of trial steps: First, a trial step in any spin direction 
uniformly distributed in spin space. This step does not 
depend on the initial direction of the spin. It samples 
the whole phase space efficiently and guarantees ergodic- 
ity. Second, a small step within a limited circular region 
around the initial spin direction. This step can efficiently 
simulate the coherent rotation. Third, a reflection of the 
spin. This step guarantees that, in the limit of large 
anisotropy, our algorithm crosses over to an efficient sim- 
ulation of an Ising-like system. For each Monte Carlo 
step we use one of these different trial steps. Our algo- 
rithm then consists of a series of Monte Carlo steps using 
the different trial steps above. Altogether, our algorithm 
is ergodic, and it guarantees that all possible reversal 
mechanisms may occur in the system and can be simu- 
lated efficiently. As we tested by comparing simulations 
with different combination s of tr ial steps, and as we will 
demonstrate in the section IIIB, for a two-spins system, 
this algorithm does not artificially change our results. 

Simulations of Heisenberg systems are much more time 
consuming, than e. g. , those of Ising systems, since the 
Heisenberg system has many more degrees of freedom. 
Apart from that, to obtain results which are comparable 
to our theoretical considerations we have to perform sim- 
ulations in the limit of low temperatures, T <C T c , where 
the critical temperature T c is 1.44 J for anisotropy d = 
p3| . Here, the metastable lifetimes are long - for single 
runs up to 5 x 10 7 MCS in our simulation. Therefore, 
for the Heisenberg system we had to restrict ourselves to 
rather small system sizes, L — 4,8, 12. However, we tried 
to minimize the statistical error by performing an aver- 
age over many Monte Carlo runs (100 ... 1000). Since 
the theoretical considerations which we want to proof 
are for finite system sizes, and since, hence, the radius 
of the particle is a variable of the theory we belive that 
the rather small system sizes of the simulation are no 
disadvantage. Apart from that, for comparison we also 
performed Monte Carlo simulations of an Ising model 
where we also used larger system sizes of up to L = 28. 

The simulations were performed on an IBM-RS6000 
workstation cluster and on two Parsytec CC parallel com- 
puters (8 and 24 PPC604 nodes, respectively). 



B. Results for the Heisenberg system 

We start our simulation with an initial spin configu- 
ration where all spins are pointing up (all spins Si — 
(0, 0, 1)). The magnetic field B = (0, 0, —B) destabilizes 
the system, and after some time the magnetization of the 
system will reverse. The metastable lifetime r is defined 
by the condition M z (r) — where M z is the z-component 
of the magnetization M = (l/N) Sj. 
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FIG. 3. Metastable lifetime r vs. 1/T for a system of two 
spins. B = 0.5J. Two different anisotropies, d — 1.5/spin 
(nucleation) and d = 0.5J/spin (coherent rotation). The solid 
lines correspond to Eq. ^| with the energy barriers explained 
in the text. 

First, we tested our algorithm by simulating the sim- 
plest imaginable system that can show both coherent ro- 
tation and nucleation - a two-spin system. Here the en- 
ergy barrier is AE n — 2(J— B) for nucleation (i. e. , the 
spins are antiparallel during the reversal) and AE cr = 
2d(l — B/(2d)) 2 for coherent rotation (i. e. , the spins 
are always parallel during the reversal). For low temper- 
atures we expect a behavior following thermal activation 
as in Eq. ^ with the energy barriers above. Figure ^ 
demonstrates that in the limit of low temperatures, we 
actually obtain constant slopes for the lnr vs. 1/T data, 
and the slopes agree perfectly with the theoretical energy 
barriers above. Hence our simulation is in agreement 
with the theoretical expectations. 

Now we turn to the simulation of larger systems. 




FIG. 4. Snapshots of simulated spin configurations at the 
lifetime r. Shown is one central plane of systems of size R — 6 
spins. B = 0.7J. a) Coherent rotation (d = 0.35J/spin, 
T = 0.09J). b) Nucleation (d = 0.7J/spin, T = 0.45J). 
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Fig. U shows spin configurations of simulated systems 
of size R = 6 spins at the metastable lifetime r. For sim- 
plicity, only one central plane of the three-dimensional 
system is shown. The z-axis of the spin components is 
pointing up. For sufficient low anisotropy (Fig.^a) the 
spins rotate nearly coherently. At the metastable lifetime 
r the magnetization vector of the system points to any 
given direction in the x — y-plane. Therefore, as horizon- 
tal component of the spins we show here that component 
of the x — y-plane of the spin space that has the largest 
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T = 0.09J). b) Nucleation (d = 0.7J/spin), T = 0.45J. 

contribution. For larger anisotropy (Figj4)b) the re- 
versal is driven by nucleation. Since it is M z (t) = 0, the 
domain wall at that time is in the center of the system di- 
viding the particle into two oppositely magnetized parts 
of equal size. 

Fig. |B| shows the corresponding time dependence of 
the z-component of the magnetization, its absolute value 

M = |M|, and its planar component M± 



Ml 



Ml 



for the same simulation from which the spin configura- 
tions of Fig stem. For the case of coherent rotation 
there is a continuous growth of the planar component of 
the magnetization during the reversal, while the abso- 
lute value of the magnetization remains nearly constant 
- apart from a small dip at the lifetime r ss 19000MCS. 
For the case of nucleation the planar component of the 
magnetization is nearly constant zero except of a small 
hump at the lifetime t « 34000MCS. Here the absolute 
value of the magnetization breaks down. 

These results lead us to the following approach to char- 



acterize the reversal mechanisms numerically: we deter- 
mine the absolute value of the magnetization at the life- 
time, M(t). In order to obtain reasonable results we have 
to take an average over many runs, so that we define a 
quantity /i — [M(t)] where the square brackets denote 
an average over many Monte Carlo runs (or systems). 
This quantity should go to zero for a nucleation-driven 
reversal, and should be finite for coherent rotation. The 
maximum value of /z in the limit of low anisotropy should 
be the spontaneous magnetization. 

In order to confirm our theoretical results numerically, 
we simulated /i for different values of the anisotropy 
d = 0.2... 2 J/spin, the magnetic field B = 0.7... 1.2 J, 
and the system size R = 2 . . . 6 spins. We took an aver- 
age over 100 (R = 6) to 1000 (R = 2) runs. Fig. | shows 
the results for the anisotropy dependence of \x for differ- 
ent system sizes and the lowest temperatures that have 
been simulated (T = 0.71 . . . 0.04 J depending on d). The 
influence of the temperature on the simulations will be 
discussed later in connection with Fig. H. 
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FIG. 6. fi vs. anisotropy for different system sizes. 
B = 0.7 J, T = 0.71 . . . 0.04J depending on d. 

As expected, for small anisotropy, \x tends to a finite 
limit while with increasing anisotropy the curves con- 
verge to zero. This effect is stronger the larger the sys- 
tem size is — a behavior that appears to be analogous to 
the finite size behavior of a system undergoing a phase 
transition. The corresponding finite-size scaling ansatz 
ill is 



where the scaling function jl 



for 



(15) 
oo so that 
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FIG. 7. Scaling plot of [i vs. anisotropy for different system 
sizes. Data correspond to Fig. H 



theoretical consideration in section |l| 



in the limit of infinite system size it is fj, ~ (doc — d) 13 
for d < doo. In order to test if this scaling form can 
be applied to the data shown in Fig. ^], we present a 
corresponding scaling plot in Fig. 0. The data collapse 
rather well, using doc = B/2 = 0.35J/spin, (3/v — 0.5 ± 
0.1, and \/v = 0.9 ±0.1. 

Obviously, in the limit R — > oo /i(d) behaves like an 
order parameter at a second-order phase transition: it is 
zero for d > doc and finite for d < d !Xll following \i ~ 
{doo-df. 

The fact that for infinite system size the transition oc- 
curs at = B/2 is in agreement with our theoretical 
considerations, since for R — > oo the region of thermally 
activated coherent rotation vanishes and the crossover 
from nucleation to (nonthermal) coherent rotation oc- 
curs at d c {R -> oo) = B/2 (Eq. [[J). That is, only for 
d < B/2 will the particle rotate coherently, and /j, must 
be finite. However, to what degree the crossover from nu- 
cleation to coherent rotation for infinite system size may 
be described as a phase transition must be left for future 
research. 

In the following we will restrict ourselves on systems of 
finite size. In order to differentiate numerically between 
the two reversal mechanisms for systems of finite size, 
we use a criterion that also comes from a study of phase 
transitions: we define the inflection point of the curves 
d(B) as that value d c (B,R) where the crossover from 
nucleation to coherent rotation occurs. These points are 
shown in Fig. ||. For large B they agree very well with 
the theoretical line. For lower B the numerical values 
are slightly too small. This systematic deviation might 
be due to the fact that the theoretical energy barrier for 
nucleation is overestimated assuming a = J/spin: due to 
occurrence of extended domain walls for lower anisotropy 
the domain wall energy might be reduced. This also re- 
duces the energy barrier of the nucleation process and 
consequently, here crossover to nucleation occurs earlier. 

Apart from the numerical determination of the cross 
overline d c discussed above, we also tried to compute the 
relevant energy barriers directly. During the simulation, 
temperature plays a crucial role. The larger d is, the 
larger is the energy barrier which has to be overcome by 
thermal activation and, hence, the higher the tempera- 
ture has to be during the simulation in order to obtain re- 
sults within a given computing time. On the other hand, 
we have to simulate as low temperatures as possible to 
see the behavior that is described by our theoretical con- 
siderations. Therefore, varying the anisotropy, we have 
to adjust the temperature. Fig. || shows the tempera- 
ture dependence of the metastable lifetime for different 
anisotropies. To save computer time, instead of an aver- 
age of the lifetimes of the individual Monte Carlo runs, 
we calculated the median Q . In the limit of low temper- 
ature we expect a behavior following thermal activation 
as in Eq. H, with the energy barrier following from the 
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FIG. 8. Metastable lifetime r vs. 1/T for four different 
anisotropies. System size R — 4 spins and B — 1J. 

As Fig. H shows, this dependence (i. e. , constant slopes 
in Fig. H for low T) can hardly be observed. All curves 
have a finite curvature even for the lowest simulated tem- 
perature except of that for d — 0.4J/spin, where the 
energy barrier is zero. We conclude that we could not 
reach low enough temperatures within the simulations 
and, hence, we analyse our data in the following way: 
We take the local slope of (In r) (1/T) as temperature- 
dependent energy barrier. These energy barriers versus 
d are shown in Fig. for three different temperatures, 
and they are compared with the theoretical results. The 
magnetic field is B = 1J. Hence, for d < 0.5 J/spin - the 
nonthermal region (Eq. |lj) - it is AE = 0. In the regime 
for thermally activated coherent rotation, i. e. , between 
d < 0.5J/spin and the crossover anisotropy d c (Eq. [i"2| ), 
the energy barrier increases following Eq. |^. Above d c 
which is roughly 0.65 J/spin, here the energy barrier re- 
mains constant since in the nucleation regime it does not 
depend on d (Eq. ||). 
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FIG. 9. Comparison of the theoretical curve for the 
d-dependence of the energy barrier with the numerical data 
for different temperatures following from Fig. ^| 

Comparing our numerical results with this theoretical 
curve, we find agreement as far as the crossover from non- 
thermal reversal to thermally driven coherent rotation is 
concerned. Also, we observe that in principal the numeri- 
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cal results seem to converge with decreasing temperature. 
However, within these simulation we cannot confirm the 
theoretical curve very accurately, especially the asymp- 
totic behavior of the nucleation regime. We conclude 
that the main reason for this quantitative deviation is 
that we do not obtain the asymptotic low-temperature 
energy barriers, since these depend still on the tempera- 
ture. 

One additional reason for deviations in the coherent 
rotation regime might be that, in order to compare the 
results of our simulation with the theoretical results, 
we had to estimate the domain- wall energy density a, 
which here we simply set to a = J/spin. This esti- 
mate might be to large for a Heisenberg system which 
can develop extended domain walls with a lower domain- 
wall energy. We could try to fit a in such a way that 
it we obtain a reasonable agreement with the numerical 
data, but this would not solve the problem mentioned 
above, namely, that we are not in the asymptotic low- 
temperature regime. 

C. Results for the Ising system 

In order to confirm our theoretical results for the en- 
ergy barrier of the nucleation regime, it is much more 
straightforward to simulate an Ising system directly in- 
stead of a Heisenberg system with large anisotropy. 
Therefore, we performed a standard Monte Carlo sim- 
ulation of an Ising system defined by the Hamiltonian 

H = -J^2 SiSj -Bj^Si (16) 

(y) * 

with Si = ±1. As before, we simulated spherical parti- 
cles with radius r = L/2 on a simple cubic lattice of size 
L x L x L with open boundary conditions. For the Ising 
system, L was varied from L = 4 to L = 32. We used the 
same methods as above, performing averages over 50-100 
Monte Carlo runs depending on the system size. 

Fig. |Io| shows the resulting temperature dependence of 
the metastable lifetime. In our data for the Ising system, 
a thermal activation corresponding to Eq^ can be much 
better extracted than for the Heisenberg system. For low 
enough temperatures, straight lines can be fitted to the 
data the slope of which determine the activation barrier 
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FIG. 10. Metastable lifetime r vs. 1/T for two different 
magnetic fields. System size R — 4 spins. 

For a comparison of our numerically determined activa- 
tion barriers with Eq|] for the theoretical energy barriers 
of the nucleation process, once more we have to estimate 
the domain-wall energy density a. In an Ising system, 
the domain- wall width is reduced to one lattice constant 
and, hence, the domain-wall energy cannot be reduced by 
an extended width of the wall. But, on the other hand, 
on a cubic lattice the energy of a domain wall per spin 
depends on the direction of the domain wall with respect 
to the axis of the lattice. This is J/spin for a wall paral- 
lel to the axis but, larger for walls in diagonal directions. 
Thus, we can expect a to be little larger than J/spin, 
and in the following it is set to a = 1.2J. 
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FIG. 11. Comparison of the theoretical curve for the energy 
barrier (Eq. ^) with the numerical data for different system 
sizes and B = 0.1 . . . 0.8. x = MBR/a. 

Fig. O compares our numerical data with Eq. |^. Since 
AE n / (Wo) depends only on the variable x = MBR/a, 
the data for different system sizes collapse on one single 
curve, and the agreement of our numerical data with the 
theoretical curve is satisfactory. 

IV. CONCLUSIONS 

We investigate the magnetization reversal of a classical 
Heisenberg system for single-domain ferromagnetic parti- 
cles. Varying the anisotropy one expects different rever- 
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sal mechanisms, the extreme cases of which are coherent 
rotation and nucleation, which in the case of a single- 
droplet-nucleation is a reversal by domain-wall motion. 
We studied the crossover from switching due to nucle- 
ation for high anisotropy, high fields, and large systems 
to coherent rotation for lower anisotropy, lower fields, and 
smaller systems. By a comparison of the relevant energy 
barriers, we derive a formula which estimates where the 
crossover from the one mechanism to the other occurs. 

If we insert the material parameters for CoPt p4[ , 
a = 0.004 J/m 2 , and d = 200kJ/m 3 in Eq. || we find 
that for vanishing magnetic field a reversal by nucleation 
should occur for particles with a radius larger than 15nm. 
While at the moment we are not aware of any direct 
measurement of this crossover, experimental hints of the 
occurrence of different, particle-size-dependent reversal 
mechanism were published by Wernsdorfer et al. [Q . 

As one important result, we found that in the limit of 
large particle size the region where a thermally activated 
coherent rotation occurs vanishes. This means that for 
large particles the rotation is always either nonthermal 
or driven by nucleation, depending only on the ratio of 
the driving field to the anisotropy. Second, in the nu- 
cleation regime the energy barrier is reduced, since here 
- in contrast to the Neel-Brown theory - for vanishing 
magnetic field the energy barrier is proportional to the 
cross-sectional area of the particle rather than its volume. 

We confirmed the result above by simulations. We also 
tried to determine the relevant energy barriers numeri- 
cally. For the case of the Heisenberg model we could 
hardly reach the low-temperatures limit which one needs 
in order to see the simplest, lowest-energy reversal mech- 
anism. However, for an Ising model, i. e. , in the limit 
of infinite anisotropy we could establish our formulas for 
the energy barrier of a single-droplet nucleation process. 
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